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ABSTRACT 

Using new data from the K2 mission, we show that WASP-47, a previously known hot Jupiter 
host, also hosts two additional transiting planets: a Neptune-sized outer planet and a super-Earth 
inner companion. We measure planetary properties from the K2 light curve and detect transit timing 
variations, confirming the planetary nature of the outer planet. We performed a large number of 
numerical simulations to study the dynamical stability of the system and to find the theoretically 
expected transit timing variations (TTVs). The theoretically predicted TTVs are in good agreement 
with those observed, and we use the TTVs to determine the masses of two planets, and place a limit 
on the third. The WASP-47 planetary system is important because companion planets can both be 
inferred by TTVs and are also detected directly through transit observations. The depth of the hot 
Jupiter’s transits make ground-based TTV measurements possible, and the brightness of the host star 
makes it amenable for precise radial velocity measurements. The system serves as a Rosetta Stone 
for understanding TTVs as a planet detection technique. 

Subject headings: planets and satellites: detection — planets and satellites: dynamical evolution and 
stability — techniques: photometric 


1. INTRODUCTION 

Due to their large sizes and short orbital periods, 
hot Jupiters (roughly Jupiter- mass p l anets with peri¬ 
ods between 0.8 and 6.3 days; iSteffen et all 12012 .1 are 
among the easiest exoplanets to detect. Both the 
first exoplanet discover ed around a main sequence star 
dMavor &; Queloa I1995T) and the firs t known transiting 
exoplanet (ICharbonneau et al.l 1200(1 iHenrv et al.l 1200(31 
were hot Jupiters. Until the launch of the Kepler space 
telescope in 2009, the majority of known transiting ex¬ 
oplanets were hot Jupiters. Hot Jupiters allow for the 
determination of many p lanetary pro perties, including 
their core masses (iBatvgin et al.ll2009f) and atmospheres 
(ICharbonneau et al.l 20021 ). For these reasons, transiting 
hot Jupiters were and continue to be the subject of many 
follow-up studies dKrcidbcrg et al.lfeOldf h 
One such follow-up study is the search for additional 
planets in the system revealed by small departures from 
perfect periodicity in the hot Jupiter transit times (called 
transit timing variations or TTVs). TTVs were pre¬ 
dicted dHolman fc Murr avl 12007)1 l Agol et ahl I2005H an d 
searched for (|Steffen fc Agoll 120051 lOibson et alJ l2009h , 
but very little evidence for TTVs was found until the 
Kepler mission discovered smaller transiting planets on 
longer perio d orbits than the ho t Jupiters detected fr om 
the ground (|Holman et al.ll20lot iLissauer et al.l[20TTI) . 

The lack of transit timing variations for hot Jupiters 
implies a dearth of nearby planets in these systems. 
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While systems exist with a known hot Jupiter and a 
distant (> 200-day pe riod) companion dKnutson et al.l 
12014 lEndl et ahl I2014f) or a warm Jupiter (orbital pe¬ 
riod 6.3 - 15.8 days) and a close-in planet (for example , 
KOI 191: ISteffen et al.ll2010t ISanchis-Oicda ct al.H2014lb 
searches for close- in, companions to hot Jupiters (as in 
ISteffen et al.ll201 2.) have not yet been successful. 

This apparent scarcity supports the idea that hot 
Jupiters form beyond the ice line and migrate inwards 
via high eccentricity migration (HEM), a process which 
would destabi li ze th e orbits of short-period companions 
(|Mustill et al.ll2015l h Studies of the Rossiter-McLaughlin 
effect have als o fou n d the fingerpr ints of high eccentric¬ 
ity migration (lAlbrccht et alJl2012lh However, statistical 
work has shown that not al l hot Jupiters can form in 
this way dDawson et al.l[2015lh so some hot Jupiters may 
have close-in planets. Addi tiona ll y, HEM may not ex¬ 
clude nearby, small planets dFogg fc Nelsonl12007 ). 

In this paper, we present an a nalysis of the WAS P-47 
system (originally announced bv lHcllier et al.l[2012f) that 
was recently observed by th e Kepler Space Tele scope in 
its new K2 operating mode dHowell et al.ll20T4) . In ad¬ 
dition to the previously known hot Jupiter in a 4.16-day 
orbit, the K2 data reveal two more transiting planets: 
a super-Earth in a 19-hour orbit, and a Neptune-sized 
planet in a 9-day orbit. We process the K2 data, deter¬ 
mine the planetary properties, and measure the transit 
times of the three planets. We find that the measured 
TTVs are consistent with the theoretical TTVs expected 
from this system and measure or place limits on the plan¬ 
ets’ masses. Finally, we perform many dynamical simu¬ 
lations of the WASP-47 system to assess its stability. 

2. K2 DATA 

Kepler observed K2 Field 3 for 69 days between 14 
November 2014 and 23 January 2015. After the data 
were publicly released, one of us (HMS) identified addi¬ 
tional transits by visual inspection of the Pre-search Data 
Conditioned (PDC) light curve of WASP-47 (designated 
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Hours from Midtransit 


Figure 1. Phase-folded short cadence K2 light curve overlaid with our best-fit transit model (red curves), and binned points (purple 
circles). In the top panel (WASP-47e), the grey circles are bins of roughly 30 seconds. In the middle and bottom panels (WASP-47b and 
WASP-47 d), the grey squares are the individual K2 short cadence datapoints. 


EPIC 206103150) produced by the Kepler/ K2 pipeline. 
We confirmed the additiona l transits by analyzing; t he K 2 
pixel level data following IVanderb urg fc Johnson! (12014H . 
A Box-Least-Squares (BLS; iKovacs et al.l 1200211 pe- 
riodogram search (as implemented in iMuirhead et all 
I2015T) of the processed long cadence light curve identified 
the 4.16-day period hot Jupiter (WASP-47 b), a Nep¬ 
tune sized planet in a 9.03-day period (WASP-47 d), and 
a super-Earth in a 0.79-day period (WASP-47 e). 

Because of the previously known hot Jupiter, WASP- 
47 was observed in K2’s “short cadence” mode, which 
consists of 58.3 second integrations in addition to the 
standard 29.41 minute “long cadence” integrations. K2 
data are dominated by systematic effects caused by 
the spacecraft’s unstable pointing which must be re¬ 
moved in order to produce high quality photometry. 
We began processing t h e sho rt cadence data following 
IVanderburg fe Johnsonl (1201411 to estimate the correla¬ 
tion between K2’s pointing and the measured flux (which 
we refer to as the K2 flat field). We used the resulting 
light curve and measured flat field as starting points in 
a simultaneous fit of the three transit signals, the flat 
field, and long term p hotometric variations (following 
IVanderburg et aLll2015H. The thre e planetary transits 
were fit with iMandel fe Agoll (120021) transit models, the 
flat field was modeled with a spline in Kepler's pointing 
position with knots placed roughly every 0.25 arcseconds, 
and the long term variations were modeled with a spline 
in time with knots placed roughly every 0.75 days. We 
performed the fit using the Levenberg-Marquardt least 


squares minimization algorithm (iMarkwardtll2009i) . The 
resulting light curvcQ shows no evidence for K2 pointing 
systematics, and yielded a photometric precision of 350 
parts per million (ppm) per 1 minute exposure. For com¬ 
parison, during its original mission, Kepler also achieved 
350 ppm per 1 minute exposure on the equally bright 
(Kp = 11.7) KOI 279. 

We measured planetary and orbital properties by fit¬ 
ting the sho rt cadence trans i t light curves of all three 
planets with IMandel fc Agoll (|2002f) transit models us¬ 
ing Markov Chain Monte Carl o (MCMC) algorith m with 
affine invariant sampling (jCoodman fc Weareil201(t ) . We 
used 50 walkers and 9000 links, and confirmed conver¬ 
gence with the test of iGewekel (j!992h and a compari¬ 
son of the Gelman-Rubin statistics for each parameter. 
We fit for the q i and <72 limb darkening parameters from 
IKippingj (|2013f l. and for each planet, we fit for the or¬ 
bital period, time of transit, orbital inclination, scaled 
semimajor axis a/i?*, and Our best-fit model is 

shown in Figure[T]and our best-fit parameters are given in 
Table [H Our measured planetary paramete rs for WASP - 
47 b a re consistent with those reported in iHellicr et all 
( 1201211 . 

We also fitted for the transit times and transit shapes 
of each transit event in the short cadence light curve si¬ 
multaneously (due to the relatively short orbital periods 
sometimes causing two transits to overlap) using MCMC. 

7 The short cadence light curve is available for download at 

www.cfa.harvard.edu/~avanderb/wasp47sc.csv 
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Figure 2. Top: observed TTVs for WASP-47e. Middle: observed 
TTVs for WASP-47 b, overlaid (for visual clarity) with a teal sine 
curve with period equal to the expected 52.67-day super-period. 
Bottom: Observed TTVs for WASP-47 d, overlaid (for visual clar¬ 
ity) with an orange sine curve with the expected super-period. 
When analyzing the transit times, we did not use the sine fits, 
they are simply to guide the eye. 


Our measured transit timcs@ are shown in Figure [2] We 
find that the TTVs of WASP-47 b and WASP-47 d are de¬ 
tected at high significance. The two TTV curves are anti¬ 
correlated and show variations on a timescale of roughly 
50 days. This is consistent with the TTV super-period 
we expect for planets orbiting in this configuration, which 
w e calculate to b e Pttv = 52.67 days using Equation 7 
of lLithwick et al.l (120121 1. 

3. VALIDATION OF WASP-47 e AND WASP-47 d 

Transiting planet signals like those we find for WASP- 
47 e and WASP-47 d can be mimicked by a variety of 
astrophysical false positive scenarios. In this section, we 
argue that this is unlikely in the case of the WASP-47 
sy stem. The ho t Ju piter, WASP-47 b, was discovered 
bv lHellier et all (12012b and confirmed with radial veloc¬ 
ity (RV) follow-up, which showed no evidence for stellar 
mass companions or spectral line shape variations, and 
detected the spectroscopic orbit of the planet. In the K2 
data, we detect transit timing variations of WASP-47 b 
which are anti-correlated with the transit timing varia¬ 
tions of WASP-47 d, and which have a super-period con¬ 
sistent with what we expect if both of these objects are 
planets. The TTVs therefore confirm that WASP-47 d 
is a planet in the same system as WASP-47 b. We also 
fitted the transit durations with a power law and found 
that they followed the expected P 1 / 3 relation (when P 
is orbital period) for planets orbiting a single star. 

The light curve is not of sufficiently high quality to 
detect TTVs for the smaller WASP-47 e, so we validate 
its pl anetary s tatus statistically. We do this using vespa 
( Mortonll2015b . an implementation of the procedure de- 


Tables available at dept.astro.lsa.umich.edu/-jcbecker/ 


scribed in iMortonl (12012b . Given constraints on back¬ 
ground sources which could be the source of the tran¬ 
sits, a constraint on the depth of any secondary eclipse, 
the host star’s parameters and location in the sky, and 
the shape of the transit light curve, vespa calculates 
the probability of a given transit signal being an astro- 
physical false positive. Both visual i nspection of ar chival 
imaging and a lucky imaging search ( IWollert et al.ll2015b 
found no close companions near WASP-47, but the lucky 
imaging is not deep enough to rule out background ob¬ 
jects that coul d cause t he shallow transits of WASP-47 e. 
Following iMontet et al.l (|2015h . we define a conservative 
radius inside of which background sources could cause 
the transits. We adopt a radius of 12 arcseconds; we de¬ 
tect the transits with photometric apertures as small as 
6 arcseconds in radius and allow for the possibility that 
stars outside of the aperture could contribute flux due 
to Kepler's 6 arcsecond point spread function. We find 
that WASP-47 e has a false positive probability (FPP) 
of roug hl y 5 x 10~ 4 . We find using the expressions from 
ILissauer et all ([2012b that since this is a three-planet sys¬ 
tem, its FPP decreases to less than 10 -5 . As such, we 
consider WASP-47 e to be validated as a bona fide planet. 

4. DYNAMICAL SIMULATIONS 
4.1. Stability Analysis 

We test the dynamical stability of the WASP-47 plan¬ 
etary system with a large ensemble of numerical simula¬ 
tions. The K2 data determine the orbital periods of the 
three bodies to high precision and place constraints on 
the other orbital elements. We sample the allowed ranges 
of the orbital elements for all three planets, randomiz¬ 
ing the orbital phases of the three bodi es. We assigned 
masse s by sampling the distribution of (Wolfgan g et al.1 
(12015b for the measured planet radii. We chose eccen¬ 
tricities from a uniform distribution that extends up to 
e = 0.3. We discard systems that do not sat i sfy th e sta¬ 
bility criteria enumerated in lFabrvckv et all (120141 ). 

Given a set of 1000 such initial conditions, we numer¬ 
ically integrate t he systems usin g the Mercury6 inte¬ 
gration package (lChamberdll999b . We use a Bulirsch- 
Stoer (B-S) integrator, requiring that system energy be 
conserved to 1 part in 10 9 . We integrate the system 
for a total simulation length of 10 Myr, unless the sys¬ 
tem goes unstable on a shorter time scale due to ejec¬ 
tion of a planet, planetary collisions, or accretion of a 
planet by the central star. To perform these computa¬ 
tionally int ensive sim ula tions, we use the Open Science 
Grid (OSG:lPordes et al.ll2007b accessed through XSEDE 
(1 Towns et al.ll20lT 

The results from this numerical survey are shown in 
Figure [3] The left panel shows the fraction of systems 
remaining stable as a function of time. About 30% of the 
systems are unstable over short time scales, and almost 
90% of the systems are unstable over long time scales. 
Once the systems reach ages of ~ 10 4 yr, they tend to 
survive over the next three orders of magnitude in inte¬ 
gration time. The remaining three panels show the mass 
and initial eccentricity of the three planets, sampled from 
the distributions specified above. One clear trend is that 
low eccentricity systems tend to survive, whereas systems 
with e p > 0.05 are generally unstable. A second trend 
that emerges from this study is that stability does not 
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Figure 3. Results from an ensemble of 1000 numerical integrations testing the stability of the system. The left panel shows the fraction of 
the systems that survive as a function of time. The other panels show the starting mass and eccentricity of the three planets sampled over 
1000 trials; the blue points represent systems that are stable, whereas the red crosses depict systems that become dynamically unstable. 


depend sensitively on the planet masses (provided that 
the orbits are nearly circular). Stable systems arise over 
a wide range of planet masses, essentially the entire range 
of masses allowed given the measured planetary radii. 

WASP-47 b and WASP-47 d orbit within about 20% of 
the 2:1 mean motion resonance (MMR). For complete¬ 
ness, we carried out a series of numerical integrations 
where the system parameters varied over the allowed, 
stable range described above. In all trials considered, 
the resonance angles were found to be circulating rather 
than librating, so there is no indication that the system 
resides in MMR. 

Each of the numerical integrations considered here 
spans 10 Myr, which corresponds to nearly one billion 
orbits of the inner planet. Tidal interactions occur on 
longer timescales than this and should be considered 
in future work. In particular, the survival of the in¬ 
ner super-Earth planet over the estimated lifetime of the 
WASP-47 system could place limits on the values of the 
tidal quality parameters Q for the bodies in the system. 

4.2. Theoretical TTVs 

We performed a second ensemble of numerical simula¬ 
tions to estimate the magnitude of transit timing varia¬ 
tions in the WASP-47 system. We used initial conditions 
similar to those adopted in the previous section, but with 
starting eccentricities e < 0.1. 

We integrated each realization of the planetary system 
for 10 years using the Mercury6 B-S integrator with time- 
steps < 0.5 seconds. We extracted transit times from 
each integration for each planet, resulting in theoreti¬ 
cal TTV curves. The resulting distributions of predicted 
TTV amplitudes are shown in Figure[4] The three distri¬ 
butions have approximately the same shape and exhibit 
well defined peaks. The TTV amplitudes we measured 
in Section [2] are consistent with the distributions we pro¬ 
duced theoretically. 

4.3. Mass Measurements from the Transit Timing 
Variations 

We measure the TTVs with high enough precision that 
dynamical fits can gi ve esti mat es of t he planetary masses. 
We use TTVFAST (|Deck et al.l I2014D to generate model 
transit times f or each observed epoch for each planet, 
and use emcee dForeman-Mackev et al.ll2013l) . an MCMC 
algorithm with affine invariant sampling, to minimize 



TTV (minutes) 


Figure 4. Histograms of potential transit timing variations for 
each planet. For a large selection of (likely) dynamically stable 
initial conditions, we integrated the system forward over a ten-year 
timescale and extracted the expected TTV amplitude. 

residuals between the observed TTVs and these model 
TTVs. In these fits, we allow each planet’s mass, eccen¬ 
tricity, argument of pericenter, and time of first transit 
to float. We imposed a uniform prior on eccentricity 
between 0 and 0.06 (as required for long-term stabil¬ 
ity). We initialized the chains with rando m arguments 
of per icenter and masses drawn from the iHellier et al.1 
(|2012D mass posterior for W ASP-47 b and the distribu¬ 
tion of I Wolfgang et al.l (I2015D for WASP-47 e and WASP- 
47 d. We used 64 walkers and 20000 iterations to explore 
the parameter space, and discarded the first 2500 itera¬ 
tions as ‘burn-in’. We confirmed t hat the MCM C chains 
had converged using the test of IGewekel (|1992l ) and the 
Gelman-Rubin statistics (which were below 1.05 for every 
parameter). The best-fit model points are overlaid with 
the observed TTVs for the outer two planets in Figure 

El 

We find that we are able to measure the masses of 
WASP-47 b and WASP-47 d and place an upper limit on 
the mass of WASP-47 e. We additionally provide limits 
on the quantities e c cos uj c — e& cos ujb- These masses and 
limits are summarized in Table [D We measure a mass of 
341^g5 M ffi for WASP-47 b, which is consistent with the 
mass measured bv IHellier et al.l (120121 ) of 362 ± 16 M® at 
the l—a level. The mass of WASP-47 d is 15.2±7 M®. 
Only an upper limit can be placed on WASP-47 e of 
<22 M®. 
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Figure 5. The best-fit theoretical points (red, with error bars) 
are overlaid with the observed TTVs (grey circles) for the best fit 
system parameters given in Tabic [l] 

5. DISCUSSION 

WASP-47 is unusual: it is the first hot Jupiter discov¬ 
ered to have additional, close-in co mpanion pla nets. Us¬ 
ing the Exoplanet Orbit Database (|Han et al.ll2014D . we 
found that among the 224 systems containing a planet 
with mass greater than 80 M® and orbital period less 
than 10 days, only six contain additional planets, and 
none of them have additional planets in orbital periods 
shorter than 100 days. That the additional planets in 
the WASP-47 system are coplanar with the hot Jupiter 
and that the planets are unstable with e > 0.05 implies 
that the WASP-47 planets either migrated in a disk or 
some damping near the end of migration took place to 
bring them into their present compact architecture. 

The continued existence of the companions in this sys¬ 
tem indicates that HEM cannot serve as the sole forma¬ 
tion mechanism for hot Jupiters. HEM would likely have 
disrupted the orbits of the smaller planets. It is quite 
possible that there is more than one potential formation 
mechanism for hot Jupiters. Additionally, recent obser¬ 
vations have identified an additional Jupiter-mass planet 
in a 571-day orbit (called WASP-47c; Neveu-VanMalle et 
al. submitted) in this system, making this the first hot 
Jupiter with both close-in companions and an external 
perturber. Future dynamical work will place limits on 
the architecture of this system. 

WASP-47 is a rare system for which planet masses 
can be determined using TTVs measured from the K2 
data set. This is because (a) the planets are far enough 
away from resonance that the super-period (52.7 days) 
is shorter than the K2 observing baseline (69 days), and 
(b) the planets are massive enough that the TTVs are 
large enough to be detectable. The detection of TTVs 
was also aided by the fact that WASP-47 was observed 
in short cadence mode, which is unusual for K2. 

Finally, WASP-47 is a favorable target for future 
follow-up observations. The V-band magnitude is 11.9, 
bright enough for precision RV follow-up studies. The 
K2 light curve shows no evidence for rotational modula¬ 
tion, indicating that WASP-47 is photometrically quiet 
and should have little RV jitter. Measuring the mass 
of the two planetary companions with RVs could both 
improve the precision of the inferred masses and test 
the consistency of TTV and RV masses, between which 


there is some tensi on (c.g. KOI 94: iWeiss et all [20131 : 
iMasuda et al.1120131) . The 1.3% depth of the transits of 
WASP-47 b makes it easily detectable from the ground. 
Previous ground based searches for transit timing varia¬ 
tions of hot Jupiters have attained timing uncertainties of 
~ 20 seconds, lo wer than the measu red TTV amplitude 
for WASP-47 b (|Gibson et al.l l2009h . Follow-up transit 
observations could place additional constraints on the 
masses of the WASP-47 planets. 

6. CONCLUSIONS 

In this work we have studied the WASP-47 planetary 
system by using data from the Kepler/ K2 mission along 
with supporting theoretical calculations. Our main re¬ 
sults can be summarized as follows: 

1. In addition to the previously known hot Jupiter 
companion WASP-47 b, the system contains two 
additional planets that are observed in transit. 
The inner planet has a ultra-short period of only 
0.789597 days, and radius of 1.829 R®. The outer 
planet has a period of 9.03081 days and a radius of 
3.63 R®, comparable to Neptune. 

2. The system is dynamically stable. We have run 
1000 10 Myr numerical integrations of the system. 
The planetary system remains stable for the 10 per¬ 
cent of the simulations that start with the lowest 
orbital eccentricities. 

3. The particular planetary system architecture of 
WASP-47 results in measurable TTVs, which are 
in good agreement with the TTVs we find from 
numerical integrations of the system. We use the 
TTVs to measure the masses of WASP-47 b (con¬ 
sistent with RV measurements) and WASP-47 d. 

4. This compact set of planets in nearly circular, 
coplanar orbits demonstrates that at least a sub¬ 
set of Jupiter-size planets can migrate in close to 
their host star in a dynamically quiet manner, sug¬ 
gesting that there may be more than one migration 
mechanism for hot Jupiters. 

The WASP-47 planetary system provides a rare oppor¬ 
tunity where planets can be both inferred from TTVs and 
seen in transit. Future observations comparing the sys¬ 
tem parameters inferred from TTVs with those inferred 
from RVs will qualitatively test TTVs as a general tech¬ 
nique. 
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Table 1 

System Parameters for WASP-47 


Parameter 

Value 


68.3% Confidence 

Comment 




Interval Width 


Stellar Parameters 

Right Ascension 

22:04:48.7 




Declination 

-12:01:08 




M* [Mq\ 

1.04 

± 

0.08 

A 

R* [*©] 

1.15 

± 

0.04 

A 

Limb darkening q\ 

0.399 

± 

0.020 

B,D 

Limb darkening <72 

0.420 

± 

0.018 

B,D 

[cgs] 

4.34 

± 

0.01 

A 

[M/H] 

+0.36 

± 

0.05 

A 

Teff [K] 

5576 

± 

67 

A 

e d cos uj d — e b cos uj b 

-0.0002 

± 

0.0195 

A,B,C 

e d sino; d - e b sin uj b 

0.0039 

± 

0.0179 

A,B,C 

WASP-47 b 

Orbital Period, P [days] 

4.1591287 

± 

0.0000049 

B 

Radius Ratio, (Rp/R*) 

0.10186 

± 

0.00023 

B 

Scaled semimajor axis, a/R* 
Orbital inclination, i [deg] 

9.715 

± 

0.050 

B 

89.03 

± 

0.27 

B 

Transit impact parameter, b 

0.164 

± 

0.045 

B 

Time of Transit tt [BJD] 

2457007.932131 

± 

0.000023 

B 

TTV amplitude [min] 

0.63 

± 

0.10 

B 

Mp [M®] 

341 


+73 

-55 

A,B,C 

Rp [Re] 

12.77 

± 

0.44 

A,B 

WASP-47 e 

Orbital Period, P [days] 

0.789597 

± 

0.000013 

B 

Radius Ratio, (Rp/R*) 

0.01456 

± 

0.00024 

B 

Scaled semimajor axis, a/R* 

3.24 

± 

0.14 

B 

Orbital inclination, i [deg] 

87.0 

± 

3.1 

B 

Transit impact parameter, b 

0.17 

± 

0.15 

B 

Time of Transit tt [BJD] 

2457011.34849 

± 

0.00038 

B 

TTV amplitude [min] 

<1.2 min for any TTV period < 80 days 

B 

Mp [M®[ 

<22 


95% Confidence 

C 

Rp [Re] 

1.829 

± 

0.070 

A,B 

WASP -47 d 





Orbital Period, P [days] 
Radius Ratio, (Rp/R*) 

9.03081 

± 

0.00019 

B 

0.02886 

± 

0.00047 

B 

Scaled semimajor axis, a/R* 
Orbital inclination, i [deg] 

16.33 

± 

0.87 

B 

89.36 

± 

0.67 

B 

Transit impact parameter, b 

0.18 

± 

0.16 

B 

Time of Transit tt [BJD] 

2457006.36927 

± 

0.00044 

B 

TTV amplitude [min] 

7.3 

± 

1.9 

B 

Mp [M®[ 

15.2 


±7 

C 

Rp [Re] 

3.63 

± 

0.14 

A,B 


Note. — A: Parameters come from IMortier et al.l (|2013l ). B: Parameters come from 
analysis of the K2 light curve. C. Parameters come from dynamical fits to the ob¬ 
served transit timing variations. D: We report the magnitude of the impact parameter, 
whereas the true value could be positive or negative. 










